Itinerant Electron Ferromagnetism in the Quantum Hall Regime 



wo: 



Marcus Kasner 1,2 , J.J. Palacios 1 ' 3 and A. H. MacDonald 1 
1 Dept. of Physics, Indiana University, Swain Hall West 117, Bloomington, IN 47405, USA 
2 Institut fur Theoretische Physik, Otto-von-Guericke Universitdt, PF 4120, D-39016 Magdeburg, Germany 
3 Department of Physics and Astronomy, University of Kentucky, Lexington KY 40506-0055, USA 

G\ 

' We report on a study of the temperature and Zeeman-coupling-strength dependence of the one- 

particle Green's function of a two-dimensional (2D) electron gas at Landau level filling factor v = 1 
where the ground state is a strong ferromagnet. Our work places emphasis on the role played by 
the itinerancy of the electrons, which carry the spin magnetization and on analogies between this 
, system and conventional itinerant electron ferromagnets. We discuss the application to this system 

of the self-consistent Hartree-Fock approximation, which is analogous to the band theory description 
OO ' of metallic ferromagnetism and fails badly at finite temperatures because it does not account for 

spin-wave excitations. We go beyond this level by evaluating the one-particle Green's function 
using a self-energy, which accounts for quasiparticle spin-wave interactions. We report results for 
the temperature dependence of the spin magnetization, the nuclear spin relaxation rate, and 2D- 
2D tunneling conductances. Our calculations predict a sharp peak in the tunneling conductance 
■ at large bias voltages with strength proportional to temperature. We compare with experiment, 

where available, and with predictions based on numerical exact diagonalization and other theoretical 
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approaches. 

PACS numbers: 73.40.Hm, 75.10.Lp, 73.20.Dx, 73.20.Mf 



I. INTRODUCTION 

The physics of a two-dimensional electron system (2DES) in a magnetic field is in many respects unique. Since 
the degeneracy of the discrete Landau levels increases in proportion to the magnetic field strength, all electrons can 
be accommodated in the lowest Landau level (LLL) for sufficiently strong fields. A Landau level then behaves, in 
many respects, like a band of zero width and the system can be regarded as the extreme limit of a strongly correlated 
narrow-band electronic system. In this paper we focus on the spin magnetism of the electron system at Landau 
' level filling factor v = 1 corresponding, when the spin degree of freedom is accounted for, to the case of a half-filled 
'Landau band'. At this filling factor, it is knownEm that, when disorder can be neglected, the ground, state of the 
2DES is ferromagnetic. We are motivated to study this system because it shares many of the difficulties!^, which have 
confounded attempts to build a complete theory of metallic ferromagnetic systems, yet is free of the troubling but 
incidental consequences of a complex band structure. The finite temperature properties of metallic ferromagnets are 
more involved than those of insulating ferromagnets because of the importance of both spin and charge degrees of 
freedom, so much so that much early theory was based on misguided attempts to assign magnetic and conducting 
properties to separate classes of electrons. Despite an immense effort and many advancescl, no completely satisfactory 
theory of metallic ferromagnets exists. Progress in understanding metallic ferromagnets has been hampered in part 
I ■ by the quantitative importance of details of the electronic band-structure, which may not be accurately known or 
may be difficult to render faithfully in going beyond mean-field theories of many-body effects. The present system 
has no such difficulties. n 

Our work is also motivated by recent experimental progress. Barrett et al.u were able to measure the temperature- 
• • . dependence of the spin-polarization and nuclear-spin relaxation rates at fixed filling factors around v — 1 using nuclear 
magnetic resonance (NMR) techniques. Later, Manfra et aln extracted the spin magnetization from magneto-optical 
absorption measurements. It is our hope that critical comparison between experiment and theory will yield insights 
with wider relevance to the finite temperature properties of itinerant electron ferromagnets. It is, however, important 
to recognize that the 2DES at v = 1 is different from conventional itinerant electron ferromagnets in several important 
ways. Most importantly, two-dimensionality implies that its spontaneous spin magnetic moment will not survive at 
finite temperatures (T c = 0.) In addition, the strong magnetic field, which through its coupling to the electron's orbital 
degrees of freedom produces Landau levels, also produces a Zeeman coupling to electron's spin degree of freedom. 
For the most studied 2DES's, those formed at GaAs/AlGaAs heterojunctions, the Zeeman coupling is quite small 
compared to both Landau level separations and the characteristic interaction energy scale. As we discuss below, the 
main effect of Zeeman coupling at low temperatures is to cut-off the decrease of the magnetization due to the thermal 
excitation of very long wavelength spin- waves and mitigate consequences of the system's reduced dimensionality. 
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The recent experimental work has stimulated two different theoretical approaches, which focus on different aspects 
of the spin-magnetization physics. Read and SachdevB have compared experimental data with large N limits of 
a quantum continuum field theory model, which provides an accurate description of the long-wavelength collective 
behavior of the electronic spins. In this theory, physical properties are dependent only on the two independent ratios 
between the thermal energy, fc^T, the Zeeman coupling strength, A 2 , and the spin-stiffness energy, p s . Recently this 
work has been extended by Timm et al.Q The field-theory description is expectedH to be accurate at low temperature 
when the Zeeman coupling strength is weak. This approach achieves reasonable overall agreement between theory 
and experiment, at least at low temperatures. Our work has somewhat different motivation and follows a different 
line. We are interested in addressing the temperature dependence of the underlying electronic structure, as it changes 
in concert with the change in the spin-magnetization. Hence we fepus on the one-particle Green's function. A brief 
account of some parts of this work has been published previously^ From the Green's function we can calculate the 
electronic spectral function and hence the magnetization, the tunneling density of states, and (if vertex corrections 
are neglected) the nuclear-spin relaxation time. The approximation we use is one, which accounts for the interaction 
of quasiparticles with thermally excited spin- waves. This approximation has deficiencies. At low-temperatures and at 
low Zeeman energies the magnetizations we calculate do not appear to be in quite as good agreement with experiment 
as magnetizations from large N approximations in the field-theory calculations. At medium and high temperatures, 
we do not account systematically for temperature-dependent screening effects which are likely to be important. Some 
progress in the latter direction has recently been reported by HaussmannE2l, whose bosonized self-consistent random- 
phase approximation yields satisfying results at high temperatures, but fails at low temperatures. Progress on these 
fronts, which can be checked by comparison with experiment, may suggest routes toward more generically satisfactory 
theories of itinerant electron ferromagnets. 

Our paper is organized as follows. In Section II we briefly review established results for the ground state and 
elementary excitations of the 2DES at v = 1, which will be important for subsequent discussion. The ground 
state has all spins aligned by an arbitrarily weak Zeeman coupling. If we neglect Landau level mixing, and we 
do throughout this paper, this state has no pure charge excitations. Its elementary excitations all have a single 
reversed spin. It turns out that, in the quantum Hall regime, the Hartree-Fock approximation (HFA) is exact for 
the ground state and the time-dependent Hartree-Fock approximation is exact for its elementary excitations. The 
situation is therefore similar to that for many typical metallic ferromagnets, where there is substantial evidence that 
the ground state is well described by the Hartree-Fock-like Kohn-Sham equationsEJ of the spin-density-functional 
formalism and that its elementary excitations are well-described by the time-dependent generalization of density- 
functional theory. In Section III we discuss the application of the self-consistent Hartree-Fock approximation (SHF) 
at finite temperature. The failure of this approximation at finite temperatures is analogous to the well known failure 
of the band theory of magnetism to provide even a rough account of the ordering temperature. The approximation, on 
which the present work is based, is discussed in Section IV. We obtain an expression for the electronic self-energy by 
analytically evaluating a particle-hole ladder summation involving Green's functions of opposite spin. We emphasize 
that in our microscopic theory, it is essential to account for screening even in the low-temperature limit. Section V 
discusses results for the spin magnetization, the spin-lattice relaxation rate, and the temperature dependent 2D-2D 
tunneling I — V relation, all based on this self-energy approximation. In this section we also compare our results 
with available experiments, with other approximate theories, and with data from finite size numerical calculations 
of the magnetization and the magnetic susceptibility. We conclude that it is necessary to account for the finite 
thickness of the quantum well in comparing with experiment and that the magnetization will be overestimated at 
high temperatures by models, which do not account for electronic itinerancy. We predict the occurrence of a sharp 
peak, with strength approximately proportional to temperature, in the tunneling I — V relation at when eV is close to 
the zero-temperature spin-splitting. In Section VI we discuss some aspects of our calculation which point to difficulties 
in developing a completely satisfactory microscopic theory. Finally, we conclude in Section VII with a brief summary. 



II. STRONG FIELD LIMIT PRELIMINARIES 



In this paper we deal with a two-dimensional electron system where the electrons are subject to a constant perpen- 
dicular magnetic field of strength B = Be z . We are particularly interested in the situation at filling factor v — 1 so 
that the number of electrons in the system (N) equals the number of single-particle orbitals available in each Landau 
level; = AB/&o- A is the area of the two-dimensional system and $o = h/e is the magnetic flux quantum. One 
of the underlying assumptions in our model is that the Landau level separation is large enough that we can ignore 
fluctuations in which electrons occupy higher orbital Landau levels. This requires that the Zeeman gap and the 
interaction strength be sufficiently small in comparison to Tiuj c where uj c = eB/m* is the cyclotron frequency. The 
microscopic Hamiltonian in the Landau gauge (A = (0, Bx, 0)) is thenEl 



2 



where A = e 2 /(47re£ c ) is the interaction coupling constant, A z = |<?/xb-B| is the Zeeman energy, and p = k y = 2irl/L y , 
(I = 0,±1, . . .) are the momenta in y-direction. The magnetic length £ c = y/h/\eB\ is used as the unit of length 
below. When Landau-level mixing is neglected the kinetic energy term in the Hamiltonian is constant (Nh~ u> c /2) and 
is therefore neglected. In Eq. (|l|), N a = J2 P c p a c p-<* ^ s ^ ne number operator for an electron with spin a. We choose 
a =| as the direction parallel to the external magnetic field. .— . 
In the Landau-gauge the two-particle matrix element is given byt_l 

W(q,p-p) = ^V(k x ,q)e-^^e ik ^-^ (2) 

where V(k) is the Fourier transform of the effective 2D Coulomb interaction, which may include modifications to 
account for the finite thickness of the quantum well containing the electrons or, as we discuss below, the ad hoc 
incorporation of screening effects. The interaction vertex, shown graphically in Fig. |l|, is a function only of k y — 
93 — 92 = Qi the momentum transfer due to the interaction, and q\ — q$ = p — p', the momentum difference of the 
interacting particles. In the physically realistic case of long-range Coulomb interactions between the electrons it is 
necessary to incorporate a neutralizing positive background by setting V(k x ,k y ) oc ^ p W{q = 0,p) = 0. Note that 
the interaction term is similar to that of a one-dimensional interacting fermion model with spin-split bands of zero 
width. However, in the present case W depends not only on the transferred momentum q, but also on the momentum 
difference of the incoming or outgoing particles. This would be true even if we choose a hard core interaction in real 
space V(r) — 47rVo<5^ 2 ^(r/£ c ). For this model, the matrix element W(q,p — p ) becomes Vo^/2/ne^^ q / 2 ~(p~p ) ) £ c/ 2 . 

Unlike most interacting electron systems, screening of mutual interactions does not play a major role in the correla- 
tion physics of a quantum Hall ferromagnet, at least at low temperatures. In fact, static screening in this limit is weak 
because of the gap for charged excitations of the ground state. Nevertheless, we obtain below the somewhat surprising 
result that the electron self-energy has a weak divergence if screening is completely neglected. For that reason we 
allow for the ad hoc inclusion of screening effects in our calculations by substituting for the Fourier transform of the 
Coulomb interaction, V c (k) = 2ir£ c /k, the local static screening form V(k) = 2ir£ c /(k + k sc ). The constant value 
assigned to the screening wave vector k sc is discussed below. 

The Green's functions we calculate at a particular temperature depend, up to an overall energy scale, only on the 
ratio A z /\. At typical field strengths this ratio is small, even though A varies approximately as B 1 / 2 while A 2 is 
proportional to B. For examplea at a magnetic field of B = 7T, A 2 /A ~ 2.2K = 0.016. Thus we will be 
interested primarily in the case where the interaction term dominates over the one-particle spin-dependent term. 

We now briefly recount some known results for the ground state and low-lying excitations of the strong field 
Hamiltonian. We start by considering the case of vanishing Zoeman coupling. If the interaction is of the hard-core 
type, the ground state at filling factor v = 1 can be determined- by identifying zero energy eigenstates of the positive 
definite Hamiltonian. It turns out that these must be the product of the the Slater determinant constructed 
from all one-particle orbitals in the lowest Landau level and a many-particle spinor. The antisymmetry property 
of the many-fermion wavefunction then requires that the spinor be completely symmetric. From this it follows that 
it has total spin quantum number S = N/2. By this argument we are able to establish with some rigor that the 
ground state at v — 1 is a strong ferromagnet, in agreemenitj with the Hartree-Fock approximation discussed in the 
following section and, if we regard the degenerate Landau level as analogous to an open atomic shell, with Hund's rule 
arguments. An infinitesimal Zeeman coupling selects from this spin multiplet a ground state, in which all spins are 
aligned; for a Zeeman field in the z direction the ground state has S z = S = N/2 and the non-degenerate ground state 
wavefunction is given exactly by the only state in the v = 1 many-fermion Hilbert space, which has these quantum 
numbers: 

|*„=i>= I TTTTTTTTT....T> |* s > . (3) 



Finite size exact diagonalization calculationsEjila can be used to establish that these conclusions remain valid for the 
realistic case of Coulomb interactions. 

The simplest neutral excitations of the ferromagnetic ground state are those with a single reversed spin_.lt turns 
out that it is possible to analytically solve for the wavefunctions and eigenenergies of these excited statesE__3. The 
normalized eigenstates may be labeled by a wavevector k and in second quantized notation are given by 
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The operator relating this state to the ground state is proportional to the Fourier transform of the projection of the 
spin-lowering operator onto the lowest Landau level: \k >= e k ^/ 4 ,§ _ (k)\^> v= i > /\/~N where S~(k) is S~(k) = 
S x (k) — iS y (k) projected onto the LLL. These states appear to be similar to the single magnon states of a localized- 
spin ferromagnet, but there is an important distinction, which is most easily explained by considering a finite size 
system. The dimension of the S z = N/2 — 1 subspace is N 2 = iV 2 , since there are 7V_possible states for the minority 
spin electron and N possible states for the majority spin hole. It is possible to showtJ that for a finite size quantum 
Hall system, the number of inequivalent values of the wavevector k is also N 2 . Since there is one wavevector for each 
state and translational invariance decouples states with different wavevectors, the eigenstates can be constructed by 
symmetry arguments alone. This should be contrasted with the case of localized spin Heisenberg ferromagnets for 
which single spin-flip states can also be constructed by symmetry arguments alone, but the number of such states 
is only N . The much larger number of states in the present problem occurs because of the possibility in itinerant 
electron systems of changing the orbital occupied by an electron whose spin has been reversed. 

The nature of the single spin-flip excitations of quantum Hall ferromagnets gradually changes from having collective 
spin wave character at long waveleng th f t o having single-particle character at larger wavevectors I 18 !! 20 ! This property 
is reflected by the dispersion relatiorJlZHla for single spin-flip states: 

e sw (k)=A z + X(a(0)-a(k)), (5) 

where the quantity^ 

S® = J-0iVWe-¥eW>**>K (6) 

Eq. (^|) is easy to understand. El It represents the attractive interaction between a minority spin electron and a 
majority spin hole, quantum mechanically smeared over their respective cyclotron orbits, and separated in real space 
by £ 2 z x k. The attractive interaction contributes negatively to the excitation energy. This magnetoexcitonic picture 
of spin-flip excitations is especially appropriate when the electron-hole separation exceeds the cyclotron orbit size, 
i.e. when £ c \k\ > 1. We use this picture in Section IV to interpret our result for the interaction between electrons 
and spin- waves. The gap for creating infinitely separated electron-hole pairs, A z + Aa(0), is associated with the 
incompressibleE3 property of quantum Hall states. The property esw(k — > 0) = A z is required by spin-rotation 
invariance of the interaction Hamiltonian. For small k the excitations are collective in character and 

e SW (k) = A 2 +Airp s £ 2 c k 2 (7) 

where p s = lim^o A(a(0) — a(k)) / (4ir£ 2 k 2 ) is the spin-stiffness parameter, which appears in field theory phenomenolo- 

gies. In Fig. |^ we plot esw(k) and the above long- wavelength approximation for the case of a weakly screened Coulomb 
interaction. 

The elementary charged excitations of the | v I / i/=i) are also known exactly: 

|fc>e = 4j*»=l) 

|*}h = Pfc,Tl**=i). (8) 

These single Slater determinant states are the maximally spin-polarized states with TV = N$ ± 1 and have energies 
E = E v= \ + fJ, + £f F and E — E v —\ — fi — £f [F , respectively. Explicit expressions for the Hartree-Fock energies are 
given in the following section. Recently it has become clear that for systems with weak Zeeman energies, charged 
Skyrmiona excitations can have lower energies than these states. However, as we discuss below, Skymion states will 
have little spectral weight in the one-particle Green's function. 



III. THE HARTREE-FOCK APPROXIMATION AND SHORTCOMINGS OF THE BAND THEORY OF 

ITINERANT ELECTRON FERROMAGNETISM 

In anticipation of subsequent sections we discuss-the Hartree-Fock approximation using the lexicon of the imaginary 
time thermodynamic Green's function technique .C3c3 When only the one-particle Zeeman term is retained in Eq. (jl]) 
the thermal Green's function (GF) is given by 
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Gi°H^n) = - -@ (9) 

where v n = (2n + 1)-k/{%0) is a fermion Matsubara frequency, (3 = 1/ksT, and £^ = —<rA z /2 — fi (a = +1 for 
|, a = — 1 for |) is the single-particle, energy measured from the chemical potential. In the strong magnetic field 
limit, translational invariance impliesEZI not only that the Green's function is diagonal in the momentum labels of the 
Landau gauge states, but also that the diagonal elements are independent of momentum. This general property leads 
to thermodynamic Green's functions, which depend only on Matsubara frequency and spin. The Feynman rules for 
the evaluation of self-energy diagrams for the present problem are summarized in Appendix 

In the diagram sum we use to approximate the electronic self-energy in the following section. The propagators, 
which appear, are self-consistent Hartree-Fock (SHF) propagators rather than the bare propagators of Eq. (^). The 
Hartree-Fock propagators are obtained by self-consistently solving the Dyson equation with the lowest order self-energy 
diagram illustrated in Fig. [| This leads to an algebraic equation for the self-energy: 

sf F = e? F - £<°» = -M0)M£? F ) ■ (10) 

Since the chemical potential at v = 1 is determined by the equation rip fa) + «f(£i) = 1 ^ follows that 
(T) = —£,^ F (T) < and that the chemical potential is fixed at /i = — A5(0)/2 independent of temperature 



The SHF Green's function is 



and Zeeman energy. For ksT < Afi(0)/4 and weak Zeeman coupling it can happen that Eq. (10) has three solu- 



tions. In this case we choose the lowest value of t;^ F since this is the solution, which minimizes the grand potential 
(fi = -2ATfc B riii(2cosh(/3£ T /2)). In Fig. | we plot the Hartree-Fock eigenvalue £f F and the spin magn etization as a 
function of temperature. The difference — ^ F = 2^ F is the exchange enhanced spin splittingEjEj of the lowest 
Landau level. Its maximum value occurs at T — and is 5,(0) A + A z , dominated by the interaction term for A 2 -C A. 
In the high temperature limit (T — > oo) this gap reduces to the bare Zeeman splitting A z . 

The sharp inflection point in the spin magnetization curve for A 2 = 0.016A in Fig. ^ is a remnant of the spontaneous 
magnetization that occurs incorrectly in the SHF at low temperatures. For v = 1 and A z = it follows from Eq. ( p^| ) 
that x = (3£j* F satisfies: 

I= « A (rRi)) ■ (12) 

At high temperatures the only solution to this equation is x — so that np(£,^ F ) — np(^ F ) = 1/2 and there is no 
spin-polarization. Expanding the r. h. s. of Eq. jl2[), we see that x ^ solutions are possible when the coefficient of 
the linear term exceeds one, i.e. for T < T^ F — a(0)A/(4fc B ). Expanding ( |l2| ) up to third order we find the expected 
mean- field behavior for T near T^ F : 

1/2 



M(T) =v^fl-^ (13) 



M(T = 0) T" F \ THF ) 



M(T = T? F ,A Z ) T f3A z \ 1/3 



Similarly, at T^ F we obtain 

M{ 

M(T = 0) ~ T" F V5(0)A 

with mean field exponent S = 3. The SHF spontaneous magnetization is plotted in Fig. ^. 

The finite T ferromagnetic instability of the SHF also appears in the the RPA expression for the spin susceptibility 
per area x ^(Qj^n) — 1/^ fo^ dr e lu) ^ T ^ ^((f, t), where u> n = 2m:/{Tij3) are bosonic Matsubara frequencies and 

X- + (q,r) = < TS-(q,r)S+faq,0) > . (15) 

The overbar is intended to emphasize that the susceptibility is to be evaluated in the strong field limit where the 
Hamiltonian can be projected onto the LLL. The RPA expression can be obtained from a ladder diagram sum with 
SHF Green's functions, similar to the sum for the self-energy detailed in the next section. The result ig_3 
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1 + I{q)X- +(HF) {l 2tt£2 (i^„ - e S w{q)) 

where the single-bubble HF spin-susceptibility is 

,hf _ HF 
2^2 {ihuJn - ^F + £HF) 



X- + W(quj)= i9flB)e — ~ 1 ] (17) 

X \1,^n) n _ n0 /. t sHF,cHF\ V 1 1 J 



is the bubble with HF lines. The effective interaction appearing in Eq. ( fl6|) is defined by 7(g) = 
-2ir£ 2 c /{g(i B ) 2 ei H2 c/ 2 a(q)\. The quantity e sw (£) 



e SW (q) = A Z + A(if * - ^f)(a(0) - a(fl) (18) 

reduces to the spin wave spectrum eswio) of Eq. (||) in the T — > limit. In this approximation, the spin- wave 
bandwidth is reduced in proportion to the spin-polarization as the temperature increases. Note that this approximation 
does not capture the finite lifetime of spin-wave states, which will result from spin-wave spin-wave interactions at 
higher temperatures. This will be one of the important limitations of the theory we present in the following section. 

For A z = the static limit of the RPA susceptibility, plotted in Fig. Q diverges at the same temperature at which the 
spontaneous magnetization determined by the SHF equations Eq. (|l(j) vanishes. These results are in disagreement with 
the Mermin- Wagner theorem, which forbids continuous broken symmetries at finite temperatures in two dimensions. 
The disagreement is expected for mean-field theory. It clear that the SHF magnetizations calculated at the relatively 
small but finite values of A z , appropriate for experimentaLpstems, will be too large. In essence, the SHF calculations 
of this section are equivalent to the Stoner band theorjoEll of metallic ferromagnetism. In both cases the ground 
state is well described. (In the present situation the SHF ground state and the exact ground state are identical.) In 
both cases the magnetization is overestimated at finite temperatures primarily because of the failure to account for 
magnetization suppression due to thermally excited collective spin- wave excitations. The approximation discussed 
in the following section remedies this gross deficiency. The situation in the present two-dimensional systems with a 
small Zccman coupling is not unlike the situation in most three-dimensional metallic ferromagnets with no external 
field. For example it is generally accepted!! 2 ] that most but not all of the magnetization suppression up to the critical 
temperature in the elemental metallic ferromagnets (Fe, Ni, and Co) is due to spin-wave excitations. Of course, 
the separation of magnetization suppression into collective spin-wave effects and the particle-hole effects cannot, in 
general, be made precise. In metallic ferromagnets, and in the present two-dimensional systems, a complete theory 
valid at moderately high temperatures requires the interplay between collective spin-fluctuations and the underlying 
fermionic degrees of freedom to be accurately described. The next section reports on an adaptation of microscopic 
theories of metallic ferromagnets to the present case. 



IV. ELECTRON SPIN- WAVE SCATTERING 



Our theory is based on an expansion in terms of SHF Green's functions. The self-energy correction to these 
propagators can in principle be expressed in terms of the exact scattering vertex r^(l, 2, 3, 4) and the exact GF. 
(Here, 1 is short for {qi,ifi n ).) We use an approximation, which, as we shall demonstrate, captures much of the 
essential .-physics. It is the analogue for quantum Hall ferromagnets of the approximation discussed by Hertz and 
EdwardslS for the case of single-band Hubbard models with ground states, which are strong ferromagnets. The self- 
energy is approximated by a particle-hole ladder summation, which gives an exact description for interactions between 
a single hole in the majority band and a single electron in the minority band. However, as the density of spin-flip 
excitations increases at higher temperatures, the approximation becomes less accurate. 

The Bethe-Salpeter integral equation for the scattering vertex is expressed diagrammatically in Fig. H. Explicit 
calculation of low order ladder-diagram particle-hole vertex parts using the Feynman rules of Appendix |A|shows that 
in this approximation 

r (4) (1, 2, 3, 4) = (q + Aq/2 - q±, Aq; nu n ) (19) 

depends only on two momenta and one frequency. Here, q = (qi + <?3)/2 and Aq = 171—93 are the center of mass 
and the relative momenta, respectively, and iuo n = i{v3, n — v \,n) is the bosonic Matsubara frequency of the relative 
motion of the ingoing and outgoing particles. 

The Bethe-Salpeter integral equation for the scattering vertex is 
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T™, (q + Aq/2 - g 4 , Aq; ioj n ) = AH>(q + Ag/2 - g 4 , Aq) 

+\Xa^{i^n)tc / - q', Aq)T%{q' + Aq/2 - q A , Aq; iuj n ). (20) 



Here, we have introduced the pair propagator, X<r,a' (* w n)i defining 

X<r,a'(iu n ) = ~ ^ (t(l/„ + U n )). (21) 

It is important to note that the g-independence of the GF in the case of quantum Hall ferromagnets immensely 
simplifies the solution of the Bethe-Salpeter equation. In fact, Eq. ( p0| ) can be reduced toj-an algebraic equation 
because the second term on the RHS is a convolution integral in the center of mass coordinate!!!! We define a partial 
Fourier transformation of the scattering function: 

/oo 
V^rit(g,Ag;^) 
-oo 

/•OO 

= e ^ Aq / 2 ' q ^£ c / dqe^T^liq + Aq/2 - g 4 , Aq; iw n ), (22) 



The corresponding transformation of the Landau-gauge vertex W(q, Aq) is the particle-hole interaction, which ap- 
peared previously in Eq. (p"8|): 

/OO p fP-n' 

dqe lpq ^W(q,Aq) = / ^^V{q')e~ q ' H "l 2 e l{q '^ q+q > )l " , (23) 
J ( 27r ) 

where we have defined the two-dimensional wavevector k = (Aq,p). With these definitions we find that 

fW (k = (p,A q );iu n ) = _ ^r-g- . (24) 

(1 - X<TM>{iu n )a{k)\) 

Our self-energy approximation consists of combining this scattering vertex with the propagators. In order to avoid 
double counting the first order terms already present in the Hartree-Fock propagators, we must subtract the first 
order term from the r. h. s. of Eq. (p4|). Inverting the transform of Eq. (22) then gives the sum of second and higher 
order terms in the Bethe-Salpeter equation. The approximate self-energy for spin a is obtained by contracting the 
incoming and outgoing lines of the opposite spin into a Hartree-Fock propagator. The result is that 

S CT (MV) = -^y^ J Ga'{i{l'n + Un))^^, (iu> n ) , (25) 



where 



The Dyson equation relating the SHF and full Green's function is 

(M^n)r 1 - {Q^^n))- 1 = -t a {iv n ). (27) 

We address below the question of whether the SHF Green's function or the corrected Green's function should be used 
in the expression for the pair propagator. If we use SHF Green's functions, the particle-hole ladder diagram of order n 
is proportional to (y^ F — v^, F At low temperature this factor is close to one and larger than the combinations 
of filling factors, which result from other diagrams of the same order. E. g., the corresponding particle-particle ladder 
is proportional to (1 — v^ F — i'^, F Y n ^ 1 \ which is zero for any T when a ^ a' . This observation may be used at low 
temperatures to systematically justify the class of diagrams aie have included. At higher temperatures we expect that 
the sums of ring diagrams, which describe screening physicsEj will be among the important omissions. 

It is possible to require that the Green's functions used in the pair-propagator and self-energy expressions be 
obtained from the Dyson equation. This leads to a set of coupled equations, which can be solved numerically. For 
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lientS 

Q-\iv n ) - {Q? F {iv n ))- 1 = -t„{iu n ) 



this purpose it is convenient!^ to express the equations in the following mixed imaginary-time imaginary-frequency 
representation: 



Ea(r)=fW,(-r)^(r) 
f[%^ n ) = 2nil [ #|( 5 ® A , - a(k)X 

Xa,a'{T) = -g a {-T)g a/ {T) (28) 

As we discuss later, we have found that these self-consistent equations tend not to have stable solutions at low 
temperature. 

Our work is based mainly on the approximation in which the pair propagator and the self-energy are evaluated 
with SHF Green's functions. The pair- propagator frequency sum may then be evaluated analytically: 



xZ>(i»n) = > ( 29 ) 



leading to the following explicit expressions for the corrections T, a (iuj n ) of the majority spin and minority spin SHF 
self-energies: 

£ T (u/ n ) = A 2 if - vf* / d(—^)a 2 {k)— — — — — -^p— 30 

and 

6 r ^ \2f hf hfs r M k 2 £ 2 c ._ 2n .{n B (esw(k)) + l-nF(^ F )} 

£;(*!/„) = A>f F - i/f F ) / d(— £) a 2 (fc) — — ' . 31 

Here, ngfiswik)) is the Bose-Einstein distribution function for the spin- waves whose dispersion is specified in Eq. (|l§|). 
For T = the occupation factors in the numerators of both self-energy expressions vanish, and the SHF-result, which 
is exact in this limit, is recovered. 

These-jelectronic self-energy expressions resemble those due to virtual phonon exchange in an electron-phonon 
system.Ej The majority spin self-energy includes contributions from processes where a majority spin electron scatters 
out to a minority spin state upon absorption of a spin-wave and processes where a minority spin electron scatters 
into a majority spin state upon emission of a spin- wave. Because the spin- wave carries spin S z = — 1 there are no 
processes where a majority spin electron scatters to a minority spin state and emits a spin-wave or a minority spin 
electron scatters to a majority spin state and absorbs a spin-wave. This distinction explains both the difference 
between the majority spin self-energy in Eq. (|3fj| ) and the minority spin self-energy in Eq. and the difference 
between these self-energies and the phonon exchange self-energy where both classes of contribution appear at onccEj. 
The electron self-energy expressions arc identical to those, which would be obtained for a model where the electrons 
and spin-waves were regarded as independent fermion and boson particles with an interaction in which fcrmions are 
scattered by emitting or absorbing spin-waves. The effective electron spin-wave interaction, which can be read-off 
from the self-energy expressions, is proportional to Xa(k). Note that, unlike the case of deformation potential electron- 
phonon coupling, the matrix clement approaches a constant as k — > 0. If the long-range Coulomb interaction is not 
screened, the electron spin-wave interaction falls off only as for large k. (For ideal 2D Coulomb interactions 

a(k) = y/n/2e~ k2 ^/ 4 I (k 2 £ 2 ./4:).) The electron spin- wave interaction at large k can be understood in terms of the 
excitonic picture discussed in Section II and is proportional to the electron-electron interaction at real space separation 
k£ 2 . If the Coulomb interaction is screened a(k) will begin to fall off more quickly once k£ 2 exceeds ~ kj^. The slow 
fall off of these large momentum transfer scattering events requires us to take account of screening when we evaluate 
the self-energy expressions. 

Note that these self-energy expressions satisfy the equation-Ef (iv n ) = — !%). This is an exact identity for 
the case v = 1, which follows from the particle-hole symmetrytij of the underlying Hamiltonian. We remark that the 
self-energies given by Eq. (|30| ) and Eq. (|3l| ) have branch cuts along a finite portion of the real line. For majority 
spins the branch cut occurs along the interval ij = (£jT = £f ,fjT = £,f F — A z ). For the minority-spin the branch 

cut interval is I[ = (£^ = + A 2 , = £,f F ))- Outside of these intervals the self-energy is real on the real line. 
Because of the branch cuts, some care is required in the numerical evaluation of the self-energy expression. 
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In this paper we concentrate on physical properties which can be expressed in terms of the one-particle real-time 
Green's function. Analytically continuing the thermal Green's function self-energy expressions, Eqs. ( ^o|) and ([h]), to 
the real frequency axis {i%v n — > E + irj) gives G r a et {E) = l/(E + if] — — E a (E)). The retarded Green's function 
is completely specified by its spectral function 

A a (E) = -2ImG r a et (E). (32) 

It is Aa-(E), which we evaluate numerically, and we start by mentioning some of its general properties. If we consider 
the system of equations defining the SHF-GF (pgj ) and start the iteration from the bottom equation with a GF 
satisfying Q-\{iv n ) — —Gi(—ivn)j we end up with a new GF satisfying the same relation, i.e. the approximate system 
of equations conserves this property. The SHF-GF has this property since = — S,f F . Since G^ et {E) — —G^"(—E), 
the spectral functions for up and down spins satisfy the following relationship: 

A 1 (E)=A l (-E). (33) 

Therefore the knowledge of the spectral function for only one spin direction for the case of v = 1 is sufficient to 
determine the result for the other spin direction. 

When the self-energy is evaluated from the Hartree-Fock GF the qualitative behavior of Af(E) can be understood 
from the analytical structure of the denominator [E — £? F — Sf (E)) of the GF. The spectral function for a =| is non- 
zero along the branch cut where the retarded self-energy has a non-zero imaginary part, i.e. for fj? < E < £^ F — A z . 
Since the real part of the self-energy is monotonically decreasing outside of this interval, it vanishes for E — > ±oo, and 
is divergent for E — > — A z from above (see below), it follows that the Green's function also has simple poles, and 
the spectral function has 5- function contributions on both sides of the branch cut. The positions of these quasiparticle 
poles in the Green's function are determined by 

E-^ F = K et (E). (34) 

Fig. H illustrates a graphical solution of this equation for T = O.lA/fes- As we discuss below, the spectral weight is 
dominated by the two 5-function contributions at E~ and Ej except at elevated temperatures. 

To get a qualitative feel for the physics of the self-energy at low temperatures, it is useful to divide it into separate 
contributions from interactions with collective long wavelength spin- waves and from interactions with short wavelength 
spin-| electron, spin-t hole pairs. We arbitrarily treat spin-waves with kl c < 1 as collective and those with k£ c > 1 
as single-particle. For collective spin-waves we can approximate the electron-spin-wave interaction by a constant 
Aa(0) and the low-energy spin-wave dispersion by A z + 4irp s £ 2 k 2 . Similarly, for the particle-hole states we can, when 
screening is neglected, approximate the interaction by \/kl c and the spin-wave energy by 2£^ F — \/k£ c - We first 

concentrate on the region of energy E near where the main spectral weight resides at low temperatures. The 
collective spin-wave contribution to the self-energy near the low energy quasiparticle is 

We have assumed here that the spin-wave Bose factor is larger than spin-down electron Fermi factor. Note that this 
requires, at a minimum, that A z be less than half the single-particle energy gap £¥ F — ^ F . The extra factor of two 
in this condition occurs because of the itinerant nature of the single-particle excitations. This self-energy contribution 
will be negative and only weakly energy-dependent for E < £? F . The particle- hole contribution to the self-energy in 
the same energy range is 

s r - (E) ,-A^>- y;^2_ E . ( 36, 

This contribution to the self-energy is sharply energy dependent, reaching a maximum for E = £¥ F , and is formally 
divergent for every energy. The divergence comes from the large k contribution to the integral, i.e. from the interaction 
of electrons with widely separated spin-flip electron-hole pairs. It can be cured by introducing screening into our theory 
as discussed earlier. Screening effectively cuts off the divergent integral at k ~ \jk sc l 2 c . At low T, a Thomas- Fermi 
approximation estimate would give k sc £ c ~ Aexp(— £^ F /ksT). The screening wavevector is exponentially small 
because of the gap for charged excitations and the ultraviolet cutoff is consequently exponentially large. However, 
except at E = £^ F , the particle-hole self-energy contribution depends only logarithmically on the screening vector 
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and the declining Fermi factor will result in a small contribution at low temperatures. At E = , the particle- hole 
contribution depends linearly on k~' leading to a very narrow peak near this energy at low temperatures. 

Because of the collective contribution to the self-energy, the low-energy quasiparticle pole will be shifted to energies 
below £f and away from this peak. In the end the particle-hole contribution to the self-energy is much larger 
than it would be for a system with short range interactions. Nevertheless, provided that the collective gap A z is 
much less than half the particle-hole gap, there will be a region at low-temperatures where its contribution becomes 
unimportant. Neglecting this contribution we find the low-energy pole has residue 

z=[l+£ 2 c dkkn B (e SW (k))] 1 (37) 
Jo 

and occurs at energy 

E- =£f F - A5(0)[l-z] (38) 

The electron-spin- wave interaction strength drops out of the expression for z because of its relationship to the spin-flip 
excitation energies. As we discuss below, the loss of spectral weight due to occupied spin- waves gives a magnetization 
suppression identical to that which would be obtained from a non-interacting spin-wave model with appropriate 
ultraviolet cutoffs. 

A similar calculation can be carried out for the high energy pole in the Green's function. For E = ^f F — A z + 5 
and 6 small we can again identify approximate collective and particle-hole contributions to the self-energy: 

Ef+( £ >*(M0)4)*f (39) 



and 



oo 



X™+(E) « A I x dkk^ _ v g (40) 



In this case, screening is not necessary to make the particle-hole contribution small at low temperatures. The second 
term in the denominator of Eq. (|3^) is negligible if 8 > k B T because of the exponential cut-off in the spin-wave Bose 
occupation factors; we see below that condition is satisfied at the upper pole in the Green's function. Keeping only 
the collective contribution gives a high energy pole with residue z + = z^ 1 — 1 and energy 

Ef = t? F -A, + Xa{0)[z- 1 -l]. (41) 

For 4np s » k B T » A z , it follows from Eq. @ that 

z- 1 -l=^ln(k B T/A z ). (42) 
87rp s 

which, since 167rp s = Aa(0), guarantees that S is in the assumed energy range. The approximations made above fail 
when — 1 is large, in which case the sums of the approximate residues exceeds one. When 1 — z is small, it follows 
from these calculations that the spectral weight of the Green's function is exhausted by the two poles in the Green's 
function. As spin-waves are excited, weight is shifted from the low energy pole, which is below the chemical potential, 
to the high-energy pole, which is above the chemical potential. In the limit of zero-temperature, the spectral weight 
lies entirely in the low-energy pole and Hartree-Fock theory results are recovered. 

As we see from Fig. ^ and from the above discussion, the real part of the self-energy diverges to — oo at the lower 
limit of the branch cut when screening is neglected and to +00 at its upper limit, even if screening is included. The 
divergence at E = ^ F in Fig. ^ produces only a small feature in our numerical calculations because these results 
were obtained with k sc =/= 0. However, the divergence at the upper boundary of the branch cut is clearly visible. The 
imaginary part of the GF is nonzero throughout the interval from ^ F to £^ F — A z because spin-flip excitations exist 
at all energies between A z and A z + Aa(0). The imaginary part of the retarded Green's function is negative definite, 
as illustrated in Fig. [?]. 

The explicit expression for the spectral density is 

2tt5(E-E-) 2ttS(E-E+) 



A a (E) 



I 1 dE \E=E-< I dE \E=E+\ 

+ 0(E - CMC - m ( - 2Jm ^ e<(£;)) - . (43) 

V " ; Ka ' ((E - ^ F - Re^{E)f + (ImE r a et (E)) 2 ) V ' 
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The condition that the integral of the spectral function over frequency equals one can be used to check the accuracy 
of our calculations, m 

We have previously^ published figures illustrating the spectral densities obtained from numerical evaluations of this 
self-energy expression at several temperatures. In agreement with the preceding analysis the spectrum consists of 
low energy and a high energy delta functions separated by a band associated with the branch cut of the self-energy. 
The incoherent band contribution to the spectral weight tends to be peaked toward its high energy extremum where 
the imaginary part of the self-energy has contributions from long-wavelcngth spin-waves. The spectral weight shifts 
with increasing temperature from the low-energy Hartree-Fock pole to the high-energy pole and partially to the 
intermediate energy band. The spectral weight shift can be understood in terms of a reduction in the probability that 
majority spins will be lined up with the exchange field from the ordered moment, which fluctuates when spin-waves 
are excited at finite temperatures. In Fig. ^ we plot the fraction of the spectral weight coming from these three 
contributions as a function of temperature. Notice that the incoherent band contribution grows rather slowly with 
temperature. In Fig. we have also plotted the positions of the two poles band as a function of temperature. Initially 
the two poles are separated by the zero-temperature exchange splitting gap. In the approximation we employ, the 
splitting increases at finite temperatures because of level repulsion with the continuum states. The results shown 
in Fig. are for screened interactions between the electrons; accounting for the finite width of the quantum well or 
heterojunction reduces the splittings to approximately 2/3 of these values for typical systems and screening associated 
with Landau level mixing will cause a further reduction. 

In the next section we discuss several experiments, which can test the predictions which result from this spin- wave 
exchange approximation for the self-energy. Potentially the most telling of these will be 2D-2D tunneling experiments 
which measure the spectral functions fairly directly. We do not expect perfect agreement between the present theory 
and experiment, although qualitative agreement seems certain. We expect that comparison with tunneling experiments 
will assist future theoretical progress. Tunneling spectrum measurements in typical band ferromagnets are much less 
informative because the band width (which is zero in the present problem) is comparable to the quasiparticle band 
spin-splitting SB 



V. OBSERVABLES 
A. The spin magnetization 

The spin-magnetization is proportional to the difference of the occupation probabilities for spin-up and spin-down 
electrons, M(T) = Mo(f| — v{) where Mq = N\gfis\/2 is the ground state spin-magnetization. Expressing the 
occupation probability, v a , in terms of the spectral functions gives: 

M(T) f°° dE , w , , s , , XN 
-Yf 1 = —n F (E)(A^E) - A^E)) 

^tanhG9£7/2)A T (£7) (44) 

■oo 2lT 

For non-interacting electrons A^(E) = 2ir5(E + A z /2) and M/Mq = tanh(/3A z /4). Since interactions tend to 
favor parallel spin alignment, we expect this result to be a lower bound for M/Mq. As discussed in more detail 
later, it should become accurate in both high and low temperature limits. Note that it reflects the itinerant nature 
of the electrons which carry the spin-magnetization. For localized spin-half particles M/Mq — tanh(/3A 2 /2); the 
magnetization is smaller at high temperatures in the itinerant case because of the number of many-particle states 
increases more rapidly with the number of reversed spins. In Hartree-Fock theory, A^(E) = 2ttS(E — ) and 
M/Mq = tanh(/3^f /F (/3)/2). This prediction for the spin-magnetization is illustrated in Fig. |l0|. Because of the 
exchange enhanced spin-splitting the magnetization is much larger at -fixed temperature than in the noninteracting 
case. This result grossly overestimates the magnetization because, asE3 in the band theory of metallic magnetism, 
magnetization suppression due to thermally excited spin-waves is not accounted for. These two simple results for 
M/Mq should therefore bound the exact result. 

Now let us turn to the results obtained using our self-energy approximation which accounts for electron-spin 
wave scattering; numerical results obtained using two different values of k sc are shown in Fig. [Io| As explained in 
the previous section, even though the screening wavevector must be finite to ensure convergence of the wavevector 
integrals in our self-energy expressions, we do not expect great sensitivity to its value until it becomes comparable 
to l^ 1 . We have estimated the appropriate value for k sc at two different temperatures as described below. For 
T ~ 0.09A/fcs and T ~ O.lSA/fcg, we find that k sc = O.Oll^ 1 and k sc = 0.1Z" 1 , respectively. The screening becomes 
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weaker at low temperatures as the system approaches incompressibility. The magnetization curves of Fig. 10 were 
calculated with these two fixed screening wavevectors; if the screening wavevector were allowed to be temperature 
dependent the magnetization value would be above the long-dashed curve in Fig. |h] for temperatures below 0.09 
A/fcs in Fig. [HI should interpolate between the long-dashed short-dashed curves temperatures between 0.09A/fcs and 
0.18A/&B, and below the short-dashed curve for temperatures beyond O.lSA/fcg. The comparably weak dependence 
of the magnetization on the screening wavevector is expected, limiting this source of uncertainty in our predictions. 
Fig. |h] shows that, for A 2 = 0.016A the magnetization decreases almost linearly with T over a wide range of 
temperatures between ~ O.OlA/fcs and ~ 0.2A/fcs. Over this temperature range, the portion of the spectral weight 
at E < is dominated by the low-energy pole, and the Fermi factor evaluated at this pole is still close to 1. Under 
these conditions the temperature dependence comes nearly entirely from that of the renormalization factor at the 
low energy pole so that M/M <~ (2z — 1) ~ 1 — (ksT) ln(fcs7 1 /A z )/47rp s . This temperature dependence is identical 
to what would be obtained from a non-interacting spin-wave model in an external magnetic field B = A z /\g/j,B\- 
Ignoring the logarithmic factor, this effect is sufficient to reduce M to small values for fc^T ~ 47r/? s ~ 0.3A, and that 
is roughly what we observe in Fig. [l(]. 

The accuracy of our calculation of spin-magnetization values is most reliably judged by comparing with results 
obtained by exact diagonalization of the many-particle Hamiltonian for a small number of electrons on a sphereO In 
Fig. |ll| we present results for the ideal Coulomb interaction obtained for N — 7, and N = 9 electrons on a sphere and 
compare them with the results obtained from our self-energy approximation. We can conclude from this comparison 
that our self-energy approximation overestimates the magnetization by approximately a factor of two at intermediate 
temperatures. We discuss the physics behind this behavior at greater length below. The magnetization values 
are least accurate at temperatures between ~ A z and ~ 0.1 A where finite-size effects, discussed below, have some 
importance and cause the magnetization per particle to be underestimated by finite-size calculations. Nevertheless, 
it seems clear that our simple self-energy expression results is an overestimate of the magnetization for 0.1A < fc^T < 
0.3A. It is interesting to compare these results with essentially exact results for a S = 1/2 Heisenberg model on 
a square latticed with a nearest neighbor exchange interactions whose strength has been adjusted to reproduce the 
spin- stiffness of the quantum Hall ferromagnet. (For the Coulomb model the spin stiffness is known exactly and 
has the value p s = Xy / n/2/16n.) For the magnetization, this model is accurate when long length scale collective 
magnetization degrees of freedom have dominant importance. For fc^T smaller than ~ 0.1 A, the Heisenberg model 
magnetizations are slightly larger than those of the N = 9 exact diagonalization calculations, presumably because 
they have smaller finite-size corrections and are therefore more accurate, see Fig. O. However, the microscopic and 
Heisenberg model M(T) curves are very similar for 0.1A < fcgT < 0.2A. The localized spin model M(T) must fall 
below the exact diagonalization results for the electron model at sufficiently high temperatures, where the Heisenberg 
model magnetization is twice as large as the magnetization of the itinerant electron system, but. this apparently does 
not occur until larger temperatures are reached. This figure also shows the results obtained!] by Timm et al. by 
evaluating leading 1/N corrections to the magnetizations of continuum SU(N) and O(N) models, extending earlier 
work by Read and Sachdev.tj At moderate temperatures, the results obtained using these approximation schemes are 
numerically better than those obtained with our self energy approximation; they are also somewhat unsatisfactory, 
however, since the SU(N) scheme leads to negative magnetizations at high temperatures and the O(N) scheme fails to 
capture the low-temperature non-interacting spin-wave limit. Moreover, the continuum field theory approach cannot 
address important microscopic electronic properties like the tunneling density of states. We discuss the experimental 
values of M(T), also shown in this figure, at greater length below. 

Our exact diagonalization results for the magnetization have the largest finite-size errors when A 2 is small. In 
order to examine the limit of small A z we have calculated the finite-size spin magnetic susceptibility, Xi which is 



plotted for N = 3, N = 5, N ~ 7 and N = 9 in Fig. 13. The non-interacting result, which in the thermodynamic 
limit approaches xo = N(gps) 2 /(SksT) is plotted for comparison. For very high T the Hartree-Fock theory result, 
X -1 = 8/cs(T— A5(0)/4)/7V(5/is) 2 , is approached. However, we see from Fig. |l3|that for T — T^ F , where \ diverges 
in Hartree-Fock theory, it is in fact enhanced by a factor of only ~ 2.5 compared to the non-interacting electron result. 
The inset in this figure expands the low-T behavior by plotting ln(x) vs. ln(/csT). In the thermodynamic limit this 
quantity should approachO Airp s /kBT as T — > 0. For a finite-size system the low T limit of x is the Curie susceptibility 
associated with the spin-quantum number of the finite-system ground state: x — > {g^s) 2 So(So + l) /iksT so that ln(x) 
is a linear function of ln(T) with an offset which increases with system size. The low T breaks in the susceptibility plots 
indicate that finite-size effects become important for temperatures smaller than ~ 0.1A/£:b, roughly consistent with 
the temperature below which the microscopic exact diagonalization and Heisenberg model Monte Carlo calculations 
differ in Fig. ||. 

Comparing a 1 these results we can conclude that while our self-energy approximation removes the gross failures of 
the Hartree-Fock approximation, it still overestimates the spin- magnetization at intermediate temperatures. Evidently, 
interactions between spin-waves accelerate the decrease of the magnetization with temperature. We can obtain some 



12 



corroboration of this interpretation by examining the low energy portion of the spectrum of the electronic Hamiltonian 
for N — 13 shown in Fig. [14]. These results are for electrons on the surface of a sphere. All excited states with positive 
total L z are shown. The largest value of L z for which a state occurs specifies its total angular momentum. The 
various eigenenergies are labeled by the occupation numbers of the corresponding non- interacting spin- wave states. 
Where the occupation number exceeds one, the label of the spin- wave state is repeated so that the number of labels is 
equal to the total spin- wave occupation number of the state. Single-spin wave states occur in the S = N/2 — 1 portion 
of the spectrum. The lowest energy multiple spin-wave state is the (1, 1) state where the L = 1 spin-wave is doubly 
occupied. Since there are three degenerate spin- waves with L = 1, if there were no spin- wave interactions, the six 
(1,1) two bosons states would be degenerate and would have an energy equal to twice the energy of the L = 1 single 



spin-wave state. The open horizontal bars in Fig. 14 indicate the non-interacting spin-wave energies. We see there 
that spin- wave interactions split the six (1,1) states into a five fold degenerate L = 2 level and a single L = level. 
Generally, multiple spin wave states are reduced in energy by interactions and this effect leads to a magnetization 
which decreases more rapidly with increasing temperature than would be expected if spin-wave interactions were 



neglected. This is most apparent in Fig. 14 in the effect of interactions on the energies of the fifteen expected (2,2) 
two-boson states. When interactions are included these states are split into L — 4, L = 2 and L — levels. We see 
in Fig. ^ that the L — 4 level, which has the largest degeneracy_.is lowered in energy by spin-wave interactions. 

We now turn to a comparison between theory and experiment .oB-An important source of uncertainty is introduced 
by the dependence of the effective interaction between 2D electrons^! on the width and, to a lesser degree, the height 
of the quantum well containing the electrons^ For a quantum well of width w and infinite barrier heights, the effective 
interactionca is V e ff(k) — F(k)V(k), whereE3 

- 327r 4 (e- fct " - 1) 8tt 2 3kw 

( ' W > ~ (kw^TT 2 + k 2 W 2 )) 2 + (kw^TT 2 + k 2 W 2 )) + (47T 2 + k 2 W 2 ) ( ' 

When the system is described by phenomenological field-theory or Heisenberg models, the microscopic physics enters 
only through the spin-stiffness. In Fig. |l5| we plot the spin-stiffness, p s , at v = 1 as a function of quantum well width 
for the case-of a quantum well with infinite barrier heights. The analytic expression relating spin-stiffness and effective 
interactionEj is 
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dfcfc d K//(fc)exp(-fc 2 ^/2). (46) 



The fc 3 factor in this integrand is responsible for considerable sensitivity of p s to V e ff(k) at large k, where the finite 
thickness corrections appear. For comparison the spin-splitting gap 5(0), plotted in Fig. [Tq, has a considerably weaker 



relative dependence on well width. As shown in Fig. 15 the spin-stiffness is reduced by approximately a factor of 



two compared to the zero-width 2D layer case for quantum wells with a width w = 3£ c . This is actually close to the 
typical experimental situation where the quantum well widths are ~ 30nm and the fields are ~ 10 Tesla. We remark 
that this reduction of p s reduces the temperature interval over which collective excitations are dominant and the 
field-theory and Heisenberg phenomenological models are appropriate. In Fig. [l6]we compare experimental data with 
exact diagonalization calculations of the spin-magnetization which account for finite-well thickness. For the purpose 
of this comparison, we consider the exact diagonalization results to be essentially exact for a model which neglects 
disorder and Landau level mixing. Compared to the best fits, the experimental magnetization decreases too slowly 
at low-temperatures, and too quickly at high temperatures. It seems clear that the experimental values are too low 
at the highest temperatures where they fall below even the non-interacting system magnetizations (this happens for 
fcgT ~ 0.09 in Barrett's experiment). The weaker T dependence at low T in the experimental data, could be due 
to disorder, and in particular to weak large length scale inhomogenity, which is neglected in all theoretical models 
discussed here. 

Finally, we comment on the role of skyrmions. One of the-. interesting results of the NMR-experiments was the 
experimental evidence for skyrmions at filling factors near oneu. Our diagrammatic theory is not able to account for 
skyrmions, although we know that neutral skyrmion-antiskyrmion pair excitations exist at v = 1. Their total energy 
for A z = is only half of that of the quasielectron-quasihole pair, although the energetic advantage drops quickly with 
increasing Zeeman gap (see Fig. 1 inEJ). Thus, the existence of such excitations at the upper end of the excitation 
spectrum should not dramatically alter the thermodynamics. 

B. Nuclear spin relaxation rate 

The optically pumped NMR experiments of Barrett et al. can also be used to measure the rate of nuclear spin- 
relaxation due to coupling to the electronic spins. At present, measurements of the temperature dependence of the 
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relaxation time over broad range of temperatures at v = 1 are unavailable. The Korpnga theory of nuclear spin 
relaxation in a metaled can be generalized to electrons in a quantum well with the resulted 
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T i = ^—£r^ " » . ' ■ (47) 



In this equation ^4 is the hyperfine coupling constant, O is the unit cell volume, and </>(z) is the envelope function of 
the electronic quantum well state. The influence of interactions on the relaxatioft^ate has been studied in the limit 
where disorder is relatively strong and interactions can be treated perturbatively.c3E3 If we neglect vertex corrections, 
the local response function m Eq. can be expressed in terms of the spectral function for the one-particle 
Green's function. The result for the relaxation rate of nuclei at the center of the quantum well is 

.k B T f°° dE -dn F (E) 2 



T^ = C { B, W )^- J_ og —--^Lx'ME)ME). (48) 

Inserting GaAs parametersEl for the prefactor in Eq. [48] gives C{B,w) = 0.47(B[T]) 3 / 2 /(w;[nm]) 2 Hz. For typical fields 
and quantum well widths the prefactor corresponds to a relaxation time ~ 100s. The dimensionless integral in Eq. 
which gives the relaxation rate in this unit, is plotted as a function of temperature in Fig. [l7| Existing experimental 
data cover only the low temperature limit, and are consistent with the very long relaxation times indicated here. 
Contributions to the relaxation rate come dominantly from the continuum portion of the spectral weight near zero 
energy. In our theory this is small at both low and high temperatures. Since T^ 1 measures the low-energy spin-flip 
excitations of the system, it is amenable to an analysis based on the continuum field theory model, which has also 
been used to obtain theoretical estimates of its temperature dependence.Q It should, in principle, be possible to extract 
more information about the spectral functions, including information on its behavior far from the Fermi energy where 
a microscopic theory is necessary, from the 2D-2D tunneling studies of quantum Hall ferromagnets, which we propose 
below. 



C. Tunneling Current 

Electronic spectral functions are traditionally measured bw tunneling experiments. The measuremcntEl of spectral 
functions for 2D electron systems is enabled by techniques-^ for making separate contact between nearby quantum 
wells. In the absence of a magnetic field, this technique has made it possible to measure the qu as i pa rticle lifetime 
including its dependence on temperature due to carrier-carrier scattering. For strong magnetic fieldscjO the tunneling 
current is related to the bias voltage by 

I{V)= hj^ J oc ^ME-eV)^n F (E))X^2A a (E)A a (E~eV) (49) 

Here, t is the tunneling amplitude and A the area of the 2D system. We caution that the above formula applies when 
both 2D layers are kept at filling factor v = 1 in the presence of a bias voltage. In 2D-2D tunneling, layer densities 
change with bias potential because of the finite capacitance of the double-layer system, unless a compensating gate 
voltage is applied. This issue is especially important at v = 1 because of the sensitivity of the electron system to 
density near this filling factor. Provided that there is no density change in either layer, each spin direction contributes 
equally to the current. Measurements at fractional filling factors have demonstrated a deep, wide, and only partially 
understood minima in the spectral function near zero energy. Our calculations suggest the possibility of further 
interesting findings when experiments are performed at v = 1. In Fig. [l^ we plot the dependence of the tunneling 
current on bias voltage for three temperatures. Within our theory a 5-peak with a substantial weight proportional to 
zz + ~ (1 — (M/Mq) 2 )/A appears. This peak arises from the product of the two poles in the spectral function, and 
it will occur at a temperature dependent value of eV equal, for the idealized case of a zero width quantum well, at 
the energy difference between the upper and the lower pole positions plotted in Fig. [9] Broader and much weaker 
features result from the convolution of a delta-function peak with the continuum continuum contribution, and still 
broader and weaker features from the self convolution of the continuum contributions to the spectral functions. Only 
this last contribution contributes to the linear tunneling conductance. In Fig. |l8| we have for visualization purposes 
arbitrarily replaced the delta-function contribution by a Lorentzian of width 0.01A. 

We expect that sharp peaks do occur at voltages near the exchange splitting, despite the quantitative limitations 
of our theory discussed above. In Fig. [l9] we plot the tunneling conductance, 
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as a function of temperature. The tunneling conductance is proportional to the square of the spectral function 
averaged over energy arguments less than ~ fcgT. It has therefore a temperature dependence similar to that of the 
nuclear spin relaxation rate. It vanishes at zero temperature and remains small even for fc^T substantially in excess 
of A z since the bulk of the spectral weight is shared between the high and low-energy poles and exchange enhanced 
spin splitting causes these to be well away from the chemical potential. 



VI. ODDS AND ENDS 



A. Consistency 

A second formally exact expression can be derived which relates the electronic spectral function to the magnetization. 
In the strong field limit, K = H — fiN = + canEJ be written as 



KAT.V.u) =^^e i ^[ihv n +^]aj;<' n ) CI) 



where 77 is a positive infinitesimal and the sum is over Matsubara frequencies. Using the spectral representation of 
the Matsubara Green's function and performing a contour integral then yields 

AT r 00 7 T? 

K a {T,V^) = -f / — n F {E){E + e ] )A«{E) . (52) 
Using thermodynamic identities, the magnetization can in turn be expressed in terms of K: 

M(T,V,n) = -~ / </>'( -7T I (•") 



PJo 




We do not obtain the same result for M from this expression as from the more direct expression discussed above if 
we go beyond the SHF approximation. This ambiguity is one of several consequences-of the fact that our self-energy 
approximation is defined in terms of Hartree-Fock propagators and is not conservingo 

Some partial self-consistency can be achieved by simply replacing the SHF occupation factors, ^ F , wherever they 
appear, by occupation factors calculated from the final spectral functions. At v = 1, it is sufficient to specify the 
difference of majority and minority spin filling factors, Av. Following this procedure requires that we solve an equation 



of the form Aj/ = f(Ais), where / incorporates the entire functional dependence on the r. h. s. of Eq.(44)). At the 
same time, we modify the HF quasiparticle energy accordingly: 

£t = -i(A, + A5(0)A^) . (54) 

Together with Eq. ( (44| ) this procedure defines an implicit equation for the difference Av, i.e. for the magnetization. In 
general this equation has at least one solution with < Av < Av HF . We have determined this solution as a function 
of temperatures and have found that the solution is unique. The result is shown in Fig. [l^. The smaller magnetization 
values are in better accord with experiment. The improvement probably does reflect a partial accounting of omissions 
of the elementary electron spin-wave scattering theory. However, the abrupt decrease of the magnetization toward 
non-interacting electron values at temperature values T ~ 0.06A/fc_B is certainly unphysical. Not surprisingly, this ad 
hoc procedure does not provide satisfying results. A more elaborate attempt at a self-consistent scheme, in which the 
full set of equations (^), was solved self-consistently has been explored by HaussmannEJ and proved to be equally 
unsatisfying. 

Some hints at possible routes toward a more accurate theosji, can be found in examinations of the transverse 
susceptibility of a quantum Hall ferromagnet. It is remarkabldlS that in the LLL, this quantity can be expressed 
exactly as a geometric series of of irreducible particle-hole bubbles ioJ n )) 

X (g,Wn = 7—77^7^: — r- (55) 
1 + I(q)<$>{q,iuj n ) 
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where I(q) is independent of frequency. In the generalized random-phase approximation, <&(<f, iw n ) is approximated 
by a bubble with Hartree-Fock propagators. Generally, the irreducible particle-hole bubble can be expressed in terms 
of an irreducible vertex function 7(<f; iv n f; i(y n + ui n ), J.). 

e -q 2 l 2 j2 1 _ 

$(g,icj„) = -(.9Mb) 2 2?r ^ 2 -gZ^7(g;^n T;i(^n+w„) |)^T(^n)^l(*( Jy « (56) 

c ii/ n 

In the RPA-expression for $ we have set 7 = 1 and the GF's are approximated by the SHF-GF's. Therefore it is 
not surprising that any improvement of Q beyond the SHF-GF makes a change of 7 necessary in order to satisfy the 
Goldstone theorem condition 1 = — J(0)$(0,0). 



B. Screening 

In any electronic system, screening of electron-electron interactions plays an important role in the many-particle 
physics. The simplest approximation for the dynamically screened interaction is the RPA, which, for the present 
problem, takes the form: 

\?<i ■ \ M) 2^4 

V(k,iuj n ) = — = = — - - = =; (57) 

e(k,iuj n ) 1 + V c (fc)n°(fc,iw„) k + 2ir£ c a°{k,iw n ) 

with the polarization function approximated by 



; - 1- 



In this approximation, Landau level degeneracy is responsible for polarization functions, which vanish at non-zero 
Matsubara frequencies. At iu> n — 0, we find that 

V(k, iu n = 0) = k2t2 Vc{k) (59) 

1 + v c (k)px £ CT n F (^ F )d MS¥ F )) 

In the long wavelength limit this leads to the following expression for the temperature dependent screening wavevector 
at v = 1 : 

ksc (T) = ^n F (£f F )(l - n F {tf F )) = g—^: (60) 

The inclusion of interaction effects in the HF-energies plays an important role in the temperature dependence of k sc 
implied by this equation. At low temperatures the screening wavevector, k sc ~ exp(— /ksT)/kBT, is extremely 
small. The limited utility of this screening approximation is evidenced by the discontinuous dependence of k sc on 
temperature found when Eq. (^0|) is solved numerically. 

Additional insight into screening in the static long-wavelength limit can be obtained by using Thomas-Fermi theory 
in which 

k sc (T)=2nl c \n 2 K = 2Tr£ c \(^) ' ( 61 ) 

\OfXy T,V,N c K^l 1 J T,V,N 

where k is the compressibility, n is the particle density and v — 2itl 2 c n. If we neglect the dependence of the spectral 
function on the chemical potential /i and exploit particle- hole symmetry at v = 1, we obtain 

ksC = ^JT f^^AmnFiE^l^n^E)). (62) 

This is an implicit equation for k sc since A^{E) itself depends on ^ F , which again depends on a(0) and thus on 
k sc (see @). With this equation we can use the improved spectral density of the spin wave-theory to estimate the 
approximate magnitude of wavevector k sc at a given temperature. Rather than solving self-consistently for k sc at each 
temperature, we have fixed k sc at two values and used Eq. ([32|) to find those temperatures at which these values are 
self-consistent. As discussed in the previous section, this procedure is adequate given the relatively weak sensitivity 
of our results at low and moderate temperatures to the value used for k sc . At higher temperatures, screening is likely 
to be important since stronger screening causes a flattening of the spin wave dispersion, which leads in turn to more 
magnetization suppression, smaller exchange splittings, more mobile charges, and hence still stronger screening. 
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C. High Temperature Behavior 



The high temperature expansion of the magnetization gives us some insight into the validity of our approximations 
and stresses the importance of screening at higher temperatures. However, these results are of rather theoretical 
interest because at large temperatures excitations to higher orbital Landau levels become more and more probable 
and our restriction to the LLL becomes questionable. 

In principle, the many-body perturbation theory expansion for the thermodynamic potential provides a systematic 
order by order expansion in powers of the interaction strength over the bare Zeeman gap, or the temperature or 
combinations thereof. However, the long range of the Coulomb interaction gives rise to divergent diagrams and 
complicates issues again; for example a logarithmically divergent contribution to the magnetization appears at third 
order in interaction strength, which can be traced back to the divergent second order bubble diagram for the self- 
energy. A systematic way to circumvent this problem is to expand in terms of interactions, which are screened by 
infinite order bubble diagram partial summations. Rather than pursue this line, we have attempted to gain some 
insight from perturbative expansions bx performing the expansion up to third order in (3 for the case of hard-core 
interaction model with V(q) — inVol^ c3. 

The leading term of the high temperature expansion for the spin magnetization is identical to the same term for free 
particles, i. e. the lowest order term in the expansion of M(T)/Mq = tanh(/3A z /4) = /3A z /4 + 0(/3 2 ). As mentioned 
previously, this limit shows that the Heisenberg model, which yields a magnetization that is twice of this value, fails 
once itinerancy becomes important. Our .calculations are based on the linked-cluster coupling-constant-integration 
expansion of the thermodynamic potentialEi The leading correction, which is quadratic in (3, is determined solely by 
the exchange integral 2Vq and leads to an increase of the spin magnetization. We have carried this expansion out to 
third order and find that 

M = M ^{l + ^f- + ^{V^ - f )} + 0(/3 4 ) (63) 

Note that the Vq term also contributes positively to the magnetization. Fig. |2(] shows MT/Mq vs. 1/T at high 
temperatures. We can see from this figure that beyond leading order interaction contributions become important for 
ksT smaller than « 2.5A. For comparison the results from exact diagonalizations for fixed particle numbers are given, 
which show the right slope but exhibit finite size corrections of the expected order of magnitude^. Obviously the 
importance of correlations in this system persists to high temperatures. This property is a key in validating the use of 
models, which correctly describe only collective fluctuations in magnetization, even when fcgT is not small compared 
to underlying interaction energy scales. 



VII. SUMMARY AND CONCLUSIONS 



In this paper we report on a study of the one-particle Green's function of a quantum Hall system at filling factor 
v = 1. The ground state of the two-dimensional electron system in this case is a strong ferromagnet. Many analogies 
exist between the properties and the theoretical description of this system and conventional metallic ferromagnets. 
At v = 1 the ground state and the elementary excitations of quantum Hall ferromagnets are given exactly by time- 
dependent Hartree-Fock theory. This success is analogous to the success of band theory in describing the ground state 
and both collective and particle-hole elementary excitations of band ferromagnets. At finite temperatures, however, 
we show that Hartree-Fock theory fails qualitatively for quantum Hall ferromagnets, just as band theory fails for 
metallic ferromagnets. Our work is based on an improved approximation for the electron self-energy, which describes 
the scattering of fermionic quasiparticles off the spin-wave collective excitations composed of coherent combinations of 
spin-flip particle-hole excitations. This perturbative approximation is equivalent to ones, which have been usedEJ for 
models of itinerant electron ferromagnets at finite temperatures. Here we have the advantage that complicated band- 
structures do not confuse the comparison of theory and experiment. We find that at intermediate temperatures where 
the density of spin- wave excitations is high, although our approximation gives a huge improvement over Hartree-Fock 
theory, it still overestimates the magnetization by nearly a factor of two. We attribute this failure to the neglect of 
interactions between spin-waves in our approximation. Nevertheless, we expect that the qualitative physics predicted 
by our approximation is correct. On the basis of our calculations we predict that a sharp peak will occur in 2D-2D 
tunneling current when the bias energy eV is approximately equal to the spin-splitting and that the strength of this 
peak will be approximately proportional to temperature at low T. This peak is due to fluctuations in the direction pf 
the exchange field, which separates the energy of majority and minority spin quasiparticles. We have recently arguedE3 
that in metallic ferromagnets, this mechanism is responsible for the temperature dependence of magnetoresistance 
in ferromagnetic tunnel junctions. In that case, however, the mechanism cannot be directly verified by tunneling 
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experiments because the width of the quasiparticle bands is comparable to their exchange-spin-splitting. Verification 
of the predicted effect in quantum Hall ferromagnets, therefore, has important implications for metallic magnetic 
tunnel junctions. 
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APPENDIX A: FEYNMAN RULES 



We list here the Feynman rules for the Landau gauge in the LLL we used to evaluate self-energy diagrams. For a 
self-energy diagram of nth order with (2n — 1) internal electron lines the rules areE3: 

1. Label all directed lines with y-momentum k y , Matsubara frequency iv n and spin index a as well as each of the 
n directed interaction lines with a two-dimensional wavevector k. Conserve the y-component of the momentum 
and the spin at each vertex. 

2. Assign an unperturbed GF Q„ [iv n ) to each line as well as a factor VXfc)e _fe2 ^/ 2 e jfcx ( 9l_93 ) £ c<5fc y; q 3 _ 92 to 
each vertex. Here, q\ stems from the outgoing line of one particle and q% from the ingoing line of the other one 
(see Fig. @). 

3. Sum up over all n free internal k y of the electron lines, the independent Matsubara frequencies and spin indices. 

4. Multiply the expression by the factor (^) n A~ n (— 1) F where F is the number of closed fermionic loops. 

5. At closed fermionic loops summations of the type 

occur guaranteeing the conservation of the transferred momenta. 
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q 3 = p'+q/2 q 2 =p'- q/2 



k y =q 3 -q 2 

q=P-q/2 ^ — L — q=p+q/2 

(7 
FIG. 1. Landau gauge electron-electron interaction vertex. 
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FIG. 2. Spin-wave dispersion esw{k) at zero temperature with a Zeeman gap A z — 0.016A and a static screening wavevector 
k sc — O.Ollc 1 . Note the small value of A z compared with the spin- wave bandwidth A2(0). The dashed line shows A z + 4irp s k 2 £ 2 -. 
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FIG. 3. Proper self-energy diagram in the self-consistent HF approximation. The propagator line in this diagram must be 
determined self-consistently. This approximation leads to a frequency independent self-energy and hence to a Green's function 
whose spectral weight consists of a single S— function. The first order tadpole diagram is absent because of the introduction of 
a neutralizing background charge. 




FIG. 4. Hartree-Fock eigenenergy £f F (solid line) as a function of temperature at v = 1; — —£f F because of particle-hole 
symmetry at v — 1. The magnetization M — Mq(v^ — v{) within the SHF is depicted for A z — 0.016A (dashed curve) and 
A z — 0.0A (long dashed curve), respectively. Note the finite magnetization at low T in the latter case incorrectly indicating the 
existence of an ordered phase for T below T c — a(0)A/(4fcs). The uniform static inverse susceptibility is plotted as a dotted 
line in units of (16tt^ A)/(sf/i B ) 2 for T > T c as well for T < T c . 



q=q 4 - Aq, i(C0+V) 



q = q- Aq/2 , ('(ca+v'J 





q 4 ,N 
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FIG. 5. The self-consistent integral equation for the scattering vertex r' 4 ' in the particle-hole ladder approximation. 
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EM 

FIG. 6. The real part of the self-energy E| et (_E) and the line E — £f F as a function of E at temperature T = O.lA/fcs 
(A z — 0.016A, k sc — 0.0U~ ,w = 0.0^ c ). The quasiparticle poles occur at the two energy values E^~ and , at which the 
curves intersect. Because the self-energy diverges in opposite directions as the upper and lower boundaries of the central 
branch-cut interval, two poles exist at any temperature in this approximation. 
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FIG. 7. The imaginary part of the spin up self-energy T,^ et (E) at T = O.lX/ks- This quantity is non-zero only, i.e. within 
the central interval 7j defined in the text. Note that the imaginary part tends to a finite value on the right boundary of this 
interval. The parameters are the same as in Fig. E3. 
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FIG. 8. Temperature dependence of the partitioning of the majority spin Green's function spectral weight between low and 
high energy poles and the intermediate energy continuum. 
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FIG. 9. Positions of the low and high energy poles in the majority spin Green's function G\ Bt {E) as a function of temperature 
(A* = 0.016A,fc sc = 0.01C\w = 0.0). 
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FIG. 10. Results for M(T) for A z = 0.016A and zero width of the quantum well. The solid curve is the spin magnetization 
for free electrons with the bare Zeeman energy A z . For comparison, the dot-dashed curve is the result of the SHF (k sc = O.Ol^ 1 ) 
theory. The long-dashed and short-dashed lines show the results of the present theory with two different choices for the ad 
hoc screening vector, k sc — 0.01Z" 1 and k sc — O.li^ 1 , which correspond to temperatures 0.09A/fcs and 0.18A/fcs, respectively. 
The magnetization results for a partially selfconsistent theory, which is discussed in Section VI, are shown as dotted curve 

(k sc = o.oic 1 )- 
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FIG. 11. Results for the magnetization from finite size diagonalizations (N — 9 - long-dashed line, N = 7 - dashed line) in 
the spherical geometry for the ideal, unscreened Coulomb interaction. This is compared with the result from our perturbative 
approximation for the self-energy (solid curve) where we assumed the screening wavevector k sc = O.OIZ^ 1 . The quantum well 
has zero width and the Zeeman energy is A z = 0.016A. 
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FIG. 12. Results for the continuous field theory: the dots show results of a Monte Carlo calculation. nThe other curves 
refer to the mean field theories for SU(AT) and O(N) models including 1/iV-corrections as described in Ref.u. For comparison 
the results from the exact diagonalization for nine particles on a sphere, the experimental data of Barrett et al. and a naive 
calculation assuming non-interacting spin-waves with a quadratic dispersion relation are added. 
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FIG. 13. Numerical finite-size results for the inverse magnetic susceptibility in units of 4\/((g/j,B) 2 N . In these units the 
free electron result is \~ = 2(fcsT/A) and the Hartree-Fock result is \ l = 2(fc_gT/A — 0.3133). 
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FIG. 14. Low energy positive L z excitations for thirteen electrons on a the surface of a sphere at v 
assignments of linear spin-wave occupations corresponding to the exact eigenstates. 
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FIG. 15. Spin-stiffness p s (solid line) and spin-splitting gap a(0) (long-dashed curve) for an infinite height quantum well as 
a function of well width w assuming an ideal Coulomb interaction. These results become exact in the strong field limit where 
Landau level mixing can be neglected. The scale for the spin-splitting energy is on the right and that for the spin-stiffness is 
on the left. 
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FIG. 16. Temperature dependence of the spin magnetization for nine particles as a function of temperature for quantum 
wells of different thicknesses compared with experiment!! In this case the quantity a describes the width of a Gaussian charge 
distribution instead of the hard-wall quantum well used elsewhere. The weak temperature dependence at low temperatures in 
experiment may be due to disorder. 
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FIG. 17. Nuclear spin relaxation rate as a function of temperature in units of the field and quantum well width dependent 
prefactor C(B, w) discussed in the text. The prefactor becomes 9.8 W~ 3 Hz for a magnetic field B = 7.05T and w = 3.114 

(k 3C = o.oic 1 )- 
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FIG. 18. Tunneling current-voltage relation as a function of bias voltage at v = 1 for three temperatures T = 0.03, 0.05 and 
0.1A/fcs where the 5-peaks of the spectral functions are replaced by Lorentzians with width e = 0.01A. The parameters used 
for this calculation are: A z = 0.016A, w = 3.11^, k ac = O.Oli^J 1 . 
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FIG. 19. Tunneling conductance in the limit of vanishing voltage as a function of temperature at v = 1 for the quantum 
well width w = 3.114 and A z = 0.016A (k sc = 0.01C 1 )- 
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FIG. 20. Spin magnetization to third order in 1/T normalized by the high-temperature free particle magnetization 
(A z /(4fcsT)) for the hard-core model (Vb = y / 7r/2,A z = 0.016A, w = 0.0). The resulting deviation of the nine particle 
exact diagonalization is due to finite size effects whose diminishing influence with increasing particle number is shown in the 
inset. 
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